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The INPOPlOa version has several improvements in the fitting process, the data sets used in the fit and in the 
general features of the solution. No big change was brought in the dynamics. As a consequence of these changes, the 
extrapolation capabilities of INPOPlOa are improved compared to INPOP08, especially for the Earth, Mars, Mercury 
and Saturn orbits. As for INPOP08, INPOPlOa provides to the user, positions, velocities of the planets and the 
moon and TT-TDB chebychev polynomials at http://www.imcce.fr/inpop 

1 The Data sets 

Several data sets have been added since INPOP08. The global distribution of the data used for the INPOP fit has 
changed its balance compared to INPOP06: now, more than 56% of the planetary data used in the fit are deduced 
from tracking data of spacecrafts including range deduced from tracking, VLBI angular positions and flyby normal 
points. The statistics of the obtained postfit residuals are presented in table[TJand [9] Comments are given in section 
4. Postfit residuals are plotted in figures^ H EH [6] and [9] 

1.1 Mercury 

• Two normal points deduced from the Mariner tracking data in 1974 and 1975 have been provided by JPL 
(Folkner 2010). These points have been included in the JPL DE solutions but not distributed. They bring 
constraints only in range with an accuracy of about 150 meters. This is an improvement of a factor 5 compared 
to the previous direct radar range observations of the Earth-Mercury distances. 

• Based on planetary ephemerides refitted by the Messenger navigation team to the tracking data of the Messenger 
mission during its flybys of Mercury (NAVPE004, NAVPE005 and NAVPE006 as provided by the SPICE NAIF) 
, we extracted 3 normal points giving 3 very accurate geocentric positions of Mercury. These positions are 
given in geocentric (a, 5, range) with an accuracy of about 1.5 mas, 2.5 mas and 5 meters respectively. 

These 5 points change drastically our knowledge of the Mercury orbit. Until now, only direct radar ranging on the 
Mercury surface were avaliable with an accuracy of about 800 meters. 

1.2 Mars and Venus 

Like with INPOP08, tracking data of MEX and VEX misions are provided by ESA (Morley 2009, Morley 2010) based 
on a procedure described in Fienga et al. (2009). For extrapolation tests, we decided to keep apart the last 6 months 
of MEX and VEX data in order to test the extrapolation capabilities of INPOPlOa (see section |4~2"|) . 

1.3 Saturn 

To the Cassini normal points provided by JPL over the 2005 to 2007 period and used in INPOP08, are added VLBI 
observations of the spacecraft published by (Jones et al. 2010). These new data are spread over a complementary 
period of time: 2 points before 2005 and 8 positions obtained from 2007 to 2009.5. They provide very accurate ICRF 
geocentric (a,S) with accuracies better than 1 mas. 

1.4 Jupiter, Uranus and Neptune 

Flybys data of the 3 planets obtained during several missions (Pionneer 10 and 11, Viking 1 and 2, Ulysses and 
Cassini) are also added, provided by Folkner (2009). The mean accuracy for the Jupiter flybys in range is about 2 
kilometers and the data are spread over 25 years. The range flybys for Uranus and Neptune are unique points with 
an accuracy of about several meters. These observations improve the estimations of the geocentric distances to the 
outer planets while no observation of that type were used in INPOP06 and INPOP08 adjustements. New optical data 
obtained from 2000 to 2008 by the United States Naval Observatory Flagstaff Station using the Flagstaff Astrometric 
Scanning Transit Telescope are also added for Uranus, Neptune and Pluto. 



1.5 Pluto 



Stellar occultations are taken into account in INPOPlOa by the use of measured offsets in topocentric (a, 5). These 
offsets have been provided by Sicardy (2009). The accuracy of these observations is about 45 mas for a and 20 mas 
for 6. Such accuracy is at least a factor 2 better than the modern CCD observations. 

1.6 Moon 

Lunar Laser Ranging observations from 1969 to 2010 are taken into account, provided by several sites: Cerga (Grasse, 
France), Mc Donald, MLRS1, MLRS2 (Fort Davis, Texas, USA), Haleakala (Hawaii, USA) and Apollo (New Mexico, 
USA). Their accuracy reaches today a few centimeters. Since INPOP08, observations from Apollo and Cerga have 
been added. The weights of recent Cerga's observations have been split into 2 periods (before and after 1995), as 
for MLRS1 ones (before and after 1985). 



2 The fitting procedure 

2.1 Fit of the mass of the sun for a fixed AU 

In INPOPlOa, the value of the AU is fixed to the value given in the IERS2003 convention. The GM of the sun 
is fitted to the observations with the initial conditions of planets, the densities of the asteroids and the oblatness 
coefficient J 2 of the sun. Values of the fitted parameters are given in table [3] 



2.2 Fit of a new selection of asteroids 

The selection of asteroids modeled in INPOPlOa is based on Kuchynka et al. (2010). The perturbations of 24635 
asteroids have been taken into account using an averaging of their perturbations by a ring with fixed physical 
characteristics for most of them and individual perturbations for 161 of them (see Kuchynka et al. 2010). Based on 
a study of the correlations between the asteroids, we found 30 asteroids among the most perturbing objects highly 
correlated with each others. In order to decrease the uncertainties on the mass estimations, we fixed 15 asteroid 
masses to values well determinated by other methods and found in Baer (2010). Besides these fixed values, we fit the 
146 GMs of asteroids using a BVLS algorithm (Lawson and Hanson, 1995) with large constraints on the densities. 

3 General features 

3.1 IAU2009 Current Best Estimates (CBE) 

We used the mass of planets provided by the IAU 2009 CBE lists (Luzum 2010). 



3.2 Fit of parameters related to LLR observations 

More than 180 parameters have a significant effect on LLR computations, including initial conditions, potential 
coefficients, reflectors and stations coordinates, biases,... But because of large correlations or better estimations 
obtained by other techniques, only 59 of them have been fitted. They have been choosen so that the ratio between 
the formal error and the fitted value is less than 5%. Values of parameters and formal errors are given in tables [6] 
[7] and [8] Initial values for Moon's coefficients of potential are from LP150Q gravity model described in Konopliv 
(2001), and some of them are then fitted. Initial values for Earth's coefficients of potential are from EGM96 
(http://cddis.nasa.gov/926/egm96/egm96.htmll; J2 and J3 are then fitted. Furthermore, two biases have 
been fitted: 

• the first one on Cerga's observations between December 1996 and June 1998. It was first noticed by Chapront- 
Touze et al. (2000). Its fitted value is -0.6 ns. 

• the second one is applied on Mc Donald's observations, between October 1972 and August 1975. Its fitted 
value is 2.2 ns. 



The final version of INPOPlOa was obtained after an iterative process between planetary fit and Moon fit. 



4 Results and Applications 



The table [TJ gives the planetary postfit residuals obtained with INPOP08 and INPOPlOa. Some of the data sets 
used in the INPOPlOa were not used in the INPOP08 adjustement: this explains some differences between the two 
columns of table[T] Postfit residuals are plotted in figures[3j [4] [5] and [6] The figures [9] show LLR residuals for each 
station and table[9] gives their standard deviations. 

4.1 Asteroid masses 

In table 2 and figures 1 and 2 are presented the asteroid masses obtained with INPOPlOa compared to values found 
in the literature. The major sources of comparisons are the values obtained with DE421 (Folkner et al. 2009), 
DE423 (Konopliv et al. 2010) and Baer (2010) with realistic errorbars estimated by Kuchynka (2010). On plots 1 
and 2, the comparisons of the asteroid GMs are ranked by their impact on the Earth-Mars distances over the 1990 
to 2010 period. It then appears clearly that the estimations for the most perturbing objects are quite consistent 
when the estimations of the weak perturbing objects show bigger discrepancies. This result is consistent with the 
planetary ephemerides method for asteroid mass computation. For the asteroids inducing perturbations up to 10 
meters, the differences in GMs are usually below 1-sigma or very close to 1-sigma except for 52 Europa. For this 
asteroid, which induced a maximum of 10 meters on the Earth-Mars distances, the errorbars are very large for all 
the determinations based on planetary ephemerides: we thus conclude to a bad determination of this mass based on 
the present interval of available data. In such case, one should have fixed this mass to the value estimated by other 
methods and especially by close-encounters as the one obtained by Baer et al (2008). This should be done for the 
next INPOP version. However, the good quality of the INPOPlOa postfit and extrapolated residuals already tells us 
that such modifications would not change significantly the global behaviour of the ephemerides. 

4.2 Extrapolation tests 

Table |4]gives the residuals obtained by extrapolation of the two ephemerides INPOP08 and INPOPlOa out of their 
fit interval. These residuals are plotted in figure 8. INPOPlOa is fitted over the [1914.2:2009.7] time interval 
when INPOP08 was fitted over the interval [1914.2:2008.2]. For Venus, the extrapolation residuals are the same for 
INPOP08 and INPOPlOa. For Mars, we first compare the extrapolation residuals of INPOPlOa and of INPOP08 
over the same time interval, from 2009.7 to 2010.35. Over this period, the residuals of INPOPlOa are 9 times smaller 
than those obtained with INPOP08 (see table|4]). This is obviously induced by the prolongation of the fitting interval 
and the addition of new data sets. In order to really test the improvement of the INPOP extrapolation capabilities, 
we compare the residuals obtained 6 months after the end of the fit period, on [2008.2:2008.8] for INPOP08 and 
over [2009.8:2010.2] for INPOPlOa. Over these two intervals, the INPOPlOa residuals have a clear linear drift of 
20 meters/yr while the INPOP08 residuals show a more quadratic trend, as one can see in figure 9. This pure 
linear trend of INPOPlOa extrapolation residuals can easily be removed (top right-hand plot of figure 9). This is 
an improvement compared to the residuals obtained with INPOP08. Indeed even after the correction of the linear 
trend, a quadratic signal is still remaining in the INPOP08 extrapolation residuals. This new behaviour of INPOPlOa 
extrapolation capabilities could be induced by the changes in the fit of the asteroid masses. 

4.3 Determination of sun J 2 and PPN (3 

Table [3] gives the obtained values for the mass of the Sun, the sun J2 and the interval of sensitivity of data to 
modification of PPN (3. Based on Fienga et al. (2010), we realized several fits for different values of j3 with 
a simultaneous fitting of initial conditions of planets, mass of the sun and asteroid densities. The given interval 
corresponds to values of j3 inducing changes in the postfit residuals below 5%. As one can see in table [3] we obtain 
with INPOPlOa an improved estimation of f3 mainly due to the use in the fit of the Mariner and Messenger normal 
points of Mercury. 

5 Work in progress 

A complementary work is in progress to improve the physical interpretation of the fitted masses of asteroids. New 
estimations of supplementary advances of perihelia of planets (especially Saturn) will be computed rapidely as well 
as combined estimations of the PPN parameters j3 and 7. 
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Table 1: Statistics of the residuals obtained after the INPOPlOa fit. 
of residuals obtained with INPOP08 are given. 



For comparisons, means and standard deviations 



Planet 






INPOP08 


INPOPlOa 


Type of Data 


Nbr 


Time 










Interval 


mean ± ler 


mean ± la 


Mercury 










Direct range [m] 


462 


1965-2000 


A A | A yi A 

30 ± 842 


1 ± 866 


Mariner range [m] 


2 


1974-1975 


-1000 ± 305 


-28 ± 85 


Messenger Mercury flybys 










Messenger ra [mas] 


3 


2008-2009 


1.1 ± 0.7 


0.4 ± 1.2 


Messenger de [mas] 


3 


r\ r\ r\ c\ aaaa 

2008-2009 


2.0 ± 1.9 


1.9 ± 2.1 


Messenger range [m] 


3 


a a a a aaaa 

2008-2009 


52 ± 619 


-0.6 ± 1.9 


Venus 










Direct range [km] 


/ion 

489 


1965-2000 


0.5 ± 2.3 


0.5 ± 2.2 


VEX range [m] 


22145 


A A A /~ aa-i a 

2006-2010 


1.6 ± 4.4 


-0.2 ± 3.9 


VLBI [mas] 


22 


-t aaa aaa^ 

1990-2007 


2 ± 2 


2 ± 2.5 


Mars range [m] 










MGS 


10474 


1998-2008 


-0.9 ± 1.6 


0.5 ± 1.9 


MEX 


24262 


A A A /— AA-I A 

2006-2010 


-3.5 ± 2.0 


0.0 ± 1.7 


Path 


90 


1997 


6.8 ± 12.5 


-5.0 ± 5.0 


Vkg 


1256 


1976-1982 


A ~~7 A 1 1 A A 

-27.4 ± 19.0 


-5.7 ± 35.0 


Mars VLBI [mas] 


96 


-1 AAA A A A — 7 

1989-2007 


0.5 ± 0.5 


-0.0 ± 0.4 


Jupiter flybys 










ra [mas] 


5 


1974-2000 


48.0 ± 40.0 


6 ± 5 


de [mas] 


5 


1974-2000 


-10.0 ± 50 


-13 ± 18 


range [km] 


5 


1974-2000 


-27± 55 


-0.6 ± 1.6 


Jupiter VLBI [mas] 


24 


1996-1997 


4 ± 11 


0.2 ± 11 


Jupiter Optical 










ra [mas] 


6216 


1914-2008 


20 ± 304 


-26 ± 304 


de [mas] 


6082 


1914-2008 


-44 ± 313 


-54 ± 303 


Saturn Cassini 










r i 

ra [mas] 


31 


2004-2007 


1.5 ± 4 


0.7 ± 4 


de [mas] 


31 


2004-2007 


7.0 ± 7 


6.5 ± 7 


range [m] 


31 


2004-2007 


0.5 ± 22 


0.0 ± 17 


Saturn VLBI Cassini 










r i 

ra [mas] 


10 


AAA A AAAA 

2004-2009 


0.3 ± 0.7 


0.0 ± 0.6 


de [mas] 


10 


2004-2009 


-1.2 ± 2.0 


0.1 ± 0.4 


Saturn Optical 










r i 
ra [mas] 


7824 


-1 A -1 A AAAA 

1914-2008 


-16 ± 305 


-16 ± 305 


de [mas] 


7799 


-1A-1 A AAAA 

1914-2008 


-7 ± 276 


-9 ± 276 


Uranus flybys 










ra [mas] 


1 


1986 


-90 


-30 


de [mas] 


1 


1986 


-36 


-7 


range [km] 


1 




1139 


0.080 


Uranus Optical 










ra [mas] 


4145 


1914-2008 


-44 ± 278 


-27 ± 290 


de [mas] 


4130 


1914-2008 


-38 ± 339 


-11 ± 338 


Neptune flybys 










ra [mas] 


1 


1989 


-88 


-11 


de [mas] 


1 


1989 


-48 


-10 


range [km] 


1 




2305 


0.004 


Neptune Optical 










ra [mas] 


4340 


1914-2008 


-32 ± 282 


2 ± 281 


de [mas] 


4320 


1914-2008 


-36 ± 335 


2 ± 330 


Pluto occultation 










ra [mas] 


13 


2005-2009 


-6 ± 46 


-1 ± 47 


de [mas] 


13 


2005-2009 


16 ± 30 


-2 ± 19 


Pluto Optical 










ra [mas] 


2449 


1914-2008 


353 ± 926 


38 ± 629 


de [mas] 


2463 


1914-2008 


-22 ± 524 


17 ± 536 



Table 2: Asteroids masses found in the recent literature and compared to the values estimated in INPOPlOa. (F) is 
for fixed value. 



Asteroids 


INPOPlOa 


Baer 2010 


DE421 


DE423 


1 A 1 1 

IAU 


inl2 * * Jl /f 
1U X Mq 


11) X Mq 


1 a12 , 7l T 

1U X Mq 


inl2 * , i\ /r 
ID X Mq 


1 


/i"7r ni/" i n o /i n 

475.836 ± 2.849 


a — 7 1 — r~ n n i a — 7 1 — 1 — 

475.500 ± 4.755 


/i/~n /inn i /i/^n/in 

468.488 ± 46.849 


At — 7 nnn 1 nrn 

467.900 ± 3.250 


2 


111.394 ± 2.808 


1 n r nnn 1 1 n/~n 

106.000 ± 1.060 


-i nn n ~7 n i 1 n nno 

100.979 ± 10.098 


-ino yi yi n 1 nr~r~n 

103.440 ± 2.550 


4 


-i oo -i o ~ 7 i -i mi 

133.137 ± 1.683 


1 00 mn i n n/"/" 

133.070 ± 0.266 


-ion nor - 1 ionn/i 

132.835 ± 13.284 


-1 on mn 1 n n/~n 

130.970 ± 2.060 


7 


—} —7—70 i -i -| /in 

7.772 ± 1.142 


/~ n /" n 1 -1 nnn 

6.860 ± 1.029 


1 nnn 1 n /"nn 

5.998 ± 0.600 


1 1— n 1 1 onn 

5.530 ± 1.320 


324 


/i /~ n n i n n~7n 

4.692 ± 0.379 




a nnn 1 n /inn 

4.980 ± 0.498 


r~ a n 1 n nnn 

5.340 ± 0.990 


3 


-i-i /~ n /i i -i n 1 o 

11.604 ± 1.313 


-1 /inn 1 n mn 

13.400 ± 2.278 


-1 -1 1 — 70 1 1 1 1 — 7 

11.573 ± 1.157 


-10 -inn 1 n n -1 n 

12.100 ± 0.910 


C 




"7 no vi i 1010 

7.Uo4 ± 1.212 


c. ac(~\ 1 n ncn 

6.460 ± 0.969 


4.558 ± 0.456 


c. ~7on 1 I a a c\ 

6.730 ± 1.640 


n 

8 


/i n~7n i n £ o I 

4.072 ± 0.631 


a /r n 1 n non 

4.260 ± 0.980 


t 770 1 n 1 ~70 

1.778 ± 0.178 


n m n i n vion 

2.010 ± 0.420 


9 


5.700 (F) 


1 — — 7 n 1 -1 ^nr~ 

5.700 ± 1.425 


/1 ~ 7 1 n /1 n — 7 

4.272 ± 0.427 


nnn 1 1 nnn 

3.280 ± 1.080 


10 


A A I~ f\f\ / 1 — \ 

44.500 (F) 


a a r~nn 1 n nnn 

44.500 ± 0.890 


yi n a ~i /" 1 a n yl 

40.416 ± 4.042 


yi a n ~~7 n t — 7 — 7 /~ 

44.970 ± 7.760 


11 


-i n n /" i i nnn 

1.886 ± 1.029 


nnn 1 n nnn 

3.090 ± 0.989 


n /"on 1 n n/~n 

2.682 ± 0.268 




13 


o nnn / 1 — \ 

8.200 (F) 


nnn 1 1 r a n 

8.200 ± 1.640 


oin/i 1 noin 

3.104 ± 0.310 




14 


a -ion / 1 — \ 

4.130 (F) 


yi -1 on 1 n nni 

4.130 ± 0.991 


n /" n n 1 n n/~n 

2.622 ± 0.262 


1 nin 1 nnin 

1.910 ± 0.810 


15 


-to nrr i i /" 1 -7 

18.856 ± 1.617 


1 r— ~ 7 n n 1 nn/in 

15.700 ± 0.942 


ino/in i 1 no/i 

12.342 ± 1.234 


i/iion 1 1 j nn 

14.180 ± 1.490 


16 


-i-i n i n i r~ -\ —? a 

11.212 ± 5.174 


1 1 nnn 1 n nnn 

11.000 ± 0.990 


1 /" n n i— 1 1 /" n n 

16.825 ± 1.682 


in jin 1 a a n 

12.410 ± 3.440 


18 


1.845 (F) 




nnin i nnni 

2.012 ± 0.201 




19 


r o o n / v \ 

6.380 (F) 


r n n 1 1 nni 

6.380 ± 1.021 


/inn 1 no>in 

3.489 ± 0.349 


nnn 1 n i - on 

3.200 ± 0.530 


nn 

20 


n or A / I - \ 

2.850 (F) 


n nrn 1 n nno 

2.850 ± 0.998 


nino 1 nnin 

2.193 ± 0.219 




21 


-1 r\-7r I -1 -1 -7 n 

1.275 ± 1.170 




1 r\ a —7 i n -1 n i - 

1.047 ± 0.105 




24 


r-\ on/" i -t n n O 

2.826 ± 1.902 


r— / — 7 n 1 1 r - i~ 

5.670 ± 2.155 


no/" 1 n n a 

3.036 ± 0.304 




25 


n nnn i n nnn 

0.002 ± 0.002 




n oni 1 n non 

0.301 ± 0.030 




28 


/i /" r~ n i n nnn 

4.652 ± 0.990 




1 n/io 1 nin/i 

1.243 ± 0.124 




29 


r~ nnn / 1 - \ 

5.920 (F) 


v nnn 1 1 nn/~ 

5.920 ± 1.006 


/" on/" i n rm 

6.826 ± 0.683 


-7 yi on 1 1 /inn 

7.420 ± 1.490 


31 


o -inn / I — \ 

3.130 (F) 


1 on i 1 nnn 

3.130 ± 1.002 


n i - n n i n nrn 

8.582 ± 0.858 




41 


nnin i n o -i 

9.213 ± 2.631 




/~\ — 7 1 1 /~\ on^ 

3.971 ± 0.397 


4.240 ± 1.770 


A O 

42 


-1 or /i i i n~7n 

1.854 ± 1.070 




n /"no 1 n n/~n 

0.693 ± 0.069 




45 


n n /" n / i — \ 

2.860 (F) 


n n/~n i n ni — 7 

2.860 ± 0.057 


n nni i n nnn 

2.991 ± 0.299 




51 


n ~7n/~ i n >ini 

0.726 ± 0.421 




1 nni - i n 1 nn 

1.085 ± 0.108 




52 


/in on/i i o nr r 

42.304 ± 8.056 


nnn i n nnn 

8.290 ± 0.995 


1 n nnn 1 1 nnn 

10.202 ± 1.020 


11 1 ~7 n 1 /i nn 

11.170 ± 8.400 


60 


n a nn i n nr 

0.402 ± 0.375 




nii~n i nni/" 

0.158 ± 0.016 




63 


n nnn i -i /"nr 

2.022 ± 1.685 




n ~7/"n i n n ~ 7 ~ 7 

0.769 ± 0.077 




65 


-7 1 -TO I y| or - /" 

7.173 ± 4.256 


non i n non 

8.930 ± 0.982 


■ onn i n r~no 

5.229 ± 0.523 




78 


o nnn i n oin 

3.232 ± 2.319 




n r a r\ i c\ r\r a 

0.640 ± 0.064 




94 


-ir~0/i~7 i ii nnn 

15.847 ± 11.522 




oiin 1 nnin 

3.119 ± 0.312 




105 


1.109 (F) 




n /" /" 1 n f\r r 

0.663 ± 0.066 




107 


-innnr - \ a /~ a ~~? 

18.205 ± 4.647 


r - /~on i ni/~n 

5.630 ± 0.169 






130 


-i-i -i v o i n non 

11.152 ± 8.032 


onn 1 n 1 nn 

3.320 ± 0.199 






135 


nni~7 i nnn/i 

0.917 ± 0.884 




n i - n n 1 n nrn 

0.588 ± 0.059 




1 on 

1J9 


it oc\a i o 01 c 

5.896 ± 3.315 




1 /ii~7 1 ni/in 

1.417 ± U.142 




1 A C 

145 


2.266 (F) 




1 ino 1 nii/i 

1.138 ± 0.114 




187 


o a r\ a i -i mr 

2.484 ± 1.075 




n — 7 n 1 1 n n ~ 7 n 

0.791 ± 0.079 




192 


0.719 (F) 




n nn/" 1 n nni 

0.806 ± 0.081 
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Figure 1: Comparisons of the first 23 asteroid masses given in 10 12 solar mass estimated by different authors and 
ranked by their impact on the Earth-Mars distances over 1990 to 2010. The mass estimations of the 3 biggest 
asteroids and perturbers Ceres, Pallas and Vesta are given in table 2. "others" stands for Baer (2010). 
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Figure 2: Comparisons of the other asteroid masses given in 10 12 solar mass estimated by different authors and 
ranked by their impact on the Earth-Mars distances over 1990 to 2010. "others" stands for Baer (2010). 



Table 3: Values of paramaters obtained by the fit of INPOP08 and INPOPlOa to observations. In INPOP08, the mass 
of the sun was not fitted directly but deduced from the fitted value of AU (Fienga et al. 2009). The (F) indicates 
that the marked values are fixed in the fit. The published uncertainty (Standish 1998) of the INPOP08 GM© fixed 
value is ± 50 km 3 . s~ 2 . The equivalent value of AU deduced from the estimation of the GM© in INPOPlOa is given 
in line 6. The INPOP a for the Earth-Moon mass ratio EMRAT, the J2© and the PPN parameter /3 are the limit of 
sensitivity of the observations to the tested parameter. For the AU and the GM©, the a are deduced from the least 
squares fit. 





INPOP08 


INPOPlOa 


DE423 




± la 


± la 




EMRAT 


(81.30054 ± 0.00005) 


(81.3005700 ± 0.0000010) 


(81.3005694 ± 0.0000015) 


J2 


(1.82 ± 0.47) x 10~ 7 


(2.40 ± 0.25) x 10~ 7 




(/3 -1) x 10- 4 


(0.75 ± 1.25) 


(0.25 ± 0.75) 


(0.4 ± 2.4) 


GM© [km 3 , s- 2 ] 


132712440017.98700 (F) 


132712440055 ± 1 


132712440042 ± 10 


AU [m] 


149597870699.2 ± 0.11 


149597870691.0 (F) 




AU [m] from GM© 




149597870704.9 ± 0.3 


149597870700.0 ± 3 



Table 4: Residuals obtained by extrapolation of INPOP08 and INPOPlOa out from the fit interval. 



Planet 

Type of Data 



INPOP08 
mean ± la 



INPOPlOa 
mean ± la 



Mars range [m] 

MEX 2009.7:2010.35 



-18 ± 63 



11.3 ± 8.6 



VEX 



2009.7:2010.35 



17 ± 50 



18 ± 50 



Table 5: Masses of asteroids fitted during the fit of INPOPlOa. The o is deduced from the least squares fit. 
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Table 6: Values of dynamical parameters fitted to LLR observations. GM E mb is the sum of Earth's and Moon's 
masses, multiplied by the gravitationnal constant and is expressed in AU 3 /day 2 . C nm E are the Earth's coefficients 
of potential (without unit). t 2 ie and t 2 2e are time delays of the Earth used for tides effects and expressed in days. 
CnmM and S nm M are the Moon's coefficients of potential (without unit). (C/MR 2 ) M is the ratio between the third 
moment of inertia of the Moon, divided by its mass and the square of the mean equatorial radius (without unit). 
kiM and tm are Love number (without unit) and time delay (in day) of the Moon. Formal errors at la are given if 
the parameter is fitted. Fixed values come from Lunar gravity model LP15Q or Earth's one EGM96. 



Name 


Value 


Formal 


error (la) 


GMemb 


8.9970116036 x 10" 


10 


±7.2 


X 


10- 


-19 


C20E 


-1.0826282367 x 


10 


-3 


±1.1 


X 


10" 


-9 


C30E 


2.8862851809 x 


10 


-6 


±3.9 


X 


10" 


-8 


C40E 


1.6196215914 x 


10 


-6 










T21E 


1.2326471124 x 


10 


-2 


±9.7 


X 


10" 


-5 


1~22E 


6.9795966104 x 


10 


-3 


±7.7 


X 


10" 


-6 


C20M 


-2.0347464196 x 


10 


-4 


±2.8 


X 


10" 


-8 


C22M 


2.2400822245 x 


10 


-5 


±2.7 


X 


10" 


-9 


C30M 


-8.3743483591 x 


10 


-6 


±2.3 


X 


10" 


-8 


C31M 


3.1934585938 x 


10 


-5 


±3.7 


X 


10" 


-7 


C32M 


4.8452131770 x 


10 


-6 










C33M 


1.7262951856 x 


10 


-6 


±6.3 


X 


10" 


-9 


C40M 


9.6422863508 x 


10 


-6 










C4IM 


-5.6926874003 x 


10 


-6 










C42M 


-1.5861997683 x 


10 


-6 










C43M 


-8.1204110561 x 


10 


-8 










C44M 


-1.2739414703 x 


10 


-7 










S3IM 


3.2781092481 x 


10 


-6 


±8.8 


X 


10" 


-8 


S32M 


1.6873163963 x 


10 


-6 


±5.7 


X 


10" 


-10 


<5*33M 


-2.4855254932 x 


10 


-7 










S4IM 


1.5743934837 x 


10 


-6 










S42M 


-1.5173124037 x 


10 


-6 










S43M 


-8.0279066453 x 


10 


-7 










Sum 


8.3147478750 x 


10 


-8 










(C/MR 2 ) M 


3.9319024676 x 


10 


-1 


±4.7 


X 


10" 


-5 


&2M 


2.6248414243 x 


10 


-2 


±1.7 


X 


10" 


4 


tm 


1.8896343615 x 


10 


-1 


±1.2 


X 


10- 


-3 



Table 7: Selenocentric coordinates of reflectors, expressed in meters. 



Reflector 


Value 


Formal error (la) 




X 


1591926.238 


±1.510 


Apollo 11 


y 


690799.993 


±3.460 




z 


21002.457 


±0.306 




X 


1652725.404 


±1.140 


Apollo 14 


y 


-520893.461 


±3.600 




z 


-109731.862 


±0.318 




X 


1554675.870 


±0.276 


Apollo 15 


y 


98193.648 


±3.380 




z 


765004.763 


±0.304 




X 


1339316.595 


±1.750 


Lunakhod2 


y 


801956.645 


±2.910 




z 


756358.336 


±0.267 



Table 8: ITRF coordinates of stations at J1997.0, expressed in meters. 



Station 


Value 


Formal prror ( 1 <r 1 

1 KJ 1 1 1 1 a 1 V_ 1 1 \J 1 I J_ Ly 1 




X 


4581 fiQ2 121 


±0.003 


Cerga 


y 


556196.025 


±0.001 




z 


4389355 020 


±0.010 






— 1330781 440 


±0.012 


Mc Donald 


y 


-5328755.514 


±0.009 




z 


3235697.496 


±0.022 




X 


-1330121.099 


±0.014 


MLRS1 


y 


-5328532.296 


±0.008 




z 


3236146.576 


±0.024 




X 


-1330021.431 


±0.002 


MLRS2 


y 


-5328403.289 


±0.003 




z 


3236481.603 


±0.010 




X 


-5466000.459 


±0.011 


Haleakala (rec.) 


y 


-2404424.720 


±0.013 




z 


2242206.708 


±0.028 




X 


-1463998.834 


±0.005 


Apollo 


y 


-5166632.675 


±0.004 




z 


3435013.101 


±0.012 



Table 9: Means and standard deviations of LLR residuals for INPOPlOa solution; both are expressed in centimeters. 
Na is the total number of observations available, Nk is the number kept in fitting process, Nr is the number that 
have been rejected according to the 3er criterion (Na=Nk+Nr). 



Station 


Period 


Mean 


Std. dev. 


Na 


Nk 


Nr 


Cerga 


1987-1995 


-0.16 


6.37 


3460 


3415 


45 


Cerga 


1995-2010 


-0.03 


4.00 


4932 


4861 


71 


Cerga 


1984-1986 


8.11 


16.00 


1187 


1158 


29 


Mc Donald 


1969-1986 


0.16 


31.78 


3604 


3489 


115 


MLRS1 


1982-1985 


-7.89 


73.28 


418 


405 


13 


MLRS1 


1985-1988 


0.22 


7.30 


174 


163 


11 


MLRS2 


1988-1996 


-0.85 


4.27 


1192 


1148 


44 


MLRS2 


1996-2008 


0.56 


4.82 


2243 


1768 


475 


Haleakala 


1984-1990 


-0.44 


8.10 


770 


734 


36 


Apollo 


2006-2009 


0.08 


4.89 


642 


640 


2 




3103. B 



saoa.s 



2103.0 



sara.o 






2GDQ.D 



2QD3.D 



200 D.3 



H»OS ■.- 



S -T 





O MsssenaeHryby ranae |m| 
>: Merrier rangB[m] 


X 








y 






1B7D 1 BSD 19E3 199D 195S 2 COD 2005 2D1D 



19TC hi: 1930 19E2 



8 



8 

8 



□ red range Mercury [m] 



* m m- 



— i — i — i — i — i — i — 

1B6& 197D 1973 1 BBQ 1&3& 1990 1993 2000 



O 
n 



o 



* one-way UG5 [m| 



-r 



O one-way Palh [m] 
— I 1 1 1 T 

199B 2000 2002 2004 2D0G 



Figure 3: Postfit residuals obtained with INPOPlOa 
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Figure 4: Postfit residuals obtained with INPOPlOa 
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Figure 5: Postfit residuals obtained with INPOPlOa 
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Figure 6: Postfit residuals obtained with INPOPlOa 
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Figure 7: Extrapolation of MEX and VEX out from the fit interval. 
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Figure 9: Postfit LLR residuals with INPOPlOa for each station, expressed in centimeters. 



Asteroids INPOPlOa Baer 2010 DE421 DE423 

IAU 10 12 x M Q 10 12 x M Q 10 12 x M Q 10 12 x M 

1 475.836 ± 2.849 475.500 ± 4.755 468.488 ± 46.849 467.900 ± 3.250 

2 111.394 ± 2.808 106.000 ± 1.060 100.979 ± 10.098 103.440 ± 2.550 
4 133.137 ± 1.683 133.070 ± 0.266 132.835 ± 13.284 130.970 ± 2.060 





INPOP08 


INPOPlOa 


DE423 




± la 


± la 


± la 


EM RAT 


(81.30054 ± 0.00005) 


(81.3005700 ± 0.0000010) 


(81.3005694 ± 0.0000015) 


J2 


(1.82 ± 0.47) x 10- 7 


(2.40 ± 0.25) x 10- 7 




09 -1) x 10" 4 


(0.75 ± 1.25) 


(0.25 ± 0.75) 


(0.4 ± 2.4) 




± 5<T 


± 5a 


± la 


GM0 [km 3 , s" 2 ] 


132712440017.98700 (F) 


132712440055 ± 1 


132712440042 ± 10 


AU [m] 


149597870699.2 ± 0.55 


149597870691.0 (F) 




AU [m] from GM© 




149597870704.9 ± 0.3 


149597870700.0 ± 3 




± la 


± la 


± la 


Ceres [10 12 x M ] 


465.8 ± 4.5 


475.836 ± 2.849 


467.900 ± 3.250 


Pallas [10 12 x M Q ] 


107.6 ± 10.0 


111.394 ± 2.808 


103.440 ± 2.550 


Vesta [10 12 x M e ] 


139.2 ± 15.0 


133.137 ± 1.683 


130.970 ± 2.060 



Table 1: Statistics of the residuals obtained after the INPOPlOa fit. For comparisons, means and standard deviations 
of residuals obtained with INPOP08 are given. 



Planet 






INPOP08 


INPOPlOa 


Type of Data 


Nbr 


Time 










Interval 


mean ± la 


mean ± la 


Mercury 










Direct range [m] 


462 


1965-2000 


30 ± 842 


7 ± 866 


Mariner range [m] 


2 


1974-1975 


-1000 ± 305 


-28 ± 85 


Messenger Mercury flybys 










Messenger ra [mas] 


3 


2008-2009 


1.1 ± 0.7 


0.4 ± 1.2 


Messenger de [mas] 


3 


2008-2009 


2.0 ± 1.9 


1.9 ± 2.1 


Messenger range [m] 


3 


2008-2009 


52 ± 619 


-0.6 ± 1.9 



